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Abstract 

We review semi-inclusive charmless B decays, focussing on threshold logarithmic resumming and on univer- 
sality of QCD dynamics in radiative and semileptonic decays. 
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1 Semileptonic B decays 



The precise determination of the magnitude of the Cabibbo-Kobayashi-Maskawa (CKM) element V u b with a 
clear uncertainty remains one of the key goals of the heavy flavor physics programs, both experimentally and 
theoretically. The smallest element in the CKM mixing matrix \V u b\ plays a crucial role in the examination of 
the unitarity constraints and of the related fundamental questions. 

The charmless semileptonic B — > X U W decay channel provides a possible path for the determination of \V u b\. 
Semileptonic B decays present several advantages, such as the possibility of using the systematic framework 
provided by the Heavy Quark Effective Theory (HQET), with the additional assumption of quark- hadron 
duality. HQET is implemented through the operator product expansion (OPE) in the form of a heavy quark 
expansion ,1 . It allows to evaluate inclusive transition rates as an expansion in inverse powers of the heavy 
quark mass. Note that quark hadron duality is not derived from first principles, but it is a necessary assumptions 
for many applications of QCD. For semi-leptonic decays the property of duality is often referred as global. For 
non-leptonic decays, where the total hadronic mass is fixed, it is only the the fact that one sums over many 
hadronic states that provides an averaging (so called local duality). The success of the QCD predictions for the 
inclusive semi-hadronic r — > is T + hadrons decays widths is a strong test of global duality [2]. 

Theoretically, issues regarding the calculation of the total semileptonic partial width T(B — > X u lv) via OPE 
are well understood 0] . The semi-leptonic decay rate can be calculated as 

r = r + ^2r 2 + ^3r 3 + --- (1) 

Here the OPE is both a nonperturbative power series in 1 /m& and a perturbative expansion in a s . The lead- 
ing term Tq is the decay rate of an on-shell b quark treated within renormalization group improved perturbative 
QCD. The perturbative corrections are known to order a 2 in the strong interactions UJ. The most remarkable 
feature of Eq. is the absence of a contribution of order 1/mj,. That means that non-perturbative corrections 
are suppressed by at least two powers of the heavy quark mass. They can be expressed as matrix elements 
of higher dimension operators in HQET and parameterized by non perturbative parameters. The absence of 
contribution of order 1/mb was observed in |S] and it is often referred as Luke's theorem. 

Theoretically the total inclusive rate would allow determination of \V u b\ to better than 10%, the main 
sources of uncertainties being the uncertainty in the b quark mass and uncertainty on potential violation of 
the underlying assumption of global quark-hadron duality. However, experimentally, the much more copious 
B — > X c lv process, which has a rate about 60 times higher, does not makes feasible a measurement over the 
full phase space. 

To overcome this background, inclusive B — » X u lv measurements utilize restricted regions of phase space 
in which the B — > X C W process is kinematically highly suppressed. The background is forbidden in the regions 
of large charged lepton energy Ei > (M B — M d )/2Mb (the endpoint), low hadronic mass Mx < Mo and 
large dilepton mass q 2 > (Mb — Mo) 2 - Extraction of from such a measurement requires knowledge of 
the fraction of the total B — » X U W that lies within the utilized region of phase space, which complicates the 
theoretical issues considerably. 

Let us consider the first two kinematic regions for which the charm background is absent, that is the large 
lepton energy region, Ei > (M B — M d )/2Mb, and the small hadronic invariant mass region Mx < Mo- In 
both cases one needs to consider the following kinematic region 

E x ~ mb, m 2 x — ^QCDm b < m b . (2) 

This kinematic region has sufficient phase space for many different resonances to be produced in the final state, 
so an inclusive description of the decays is still appropriate. However, in this region, the differential rate is 
very sensitive to the details of the wave function of the b quark in the B meson. The parton level differential 
distribution at the end-point region has its own problems, as well, related to the presence of large logarithms 
which spoil the perturbative expansion. 

A third way to isolate the charmless semileptpnic signal is to use a selection based on the q 2 of the leptonic 
system. Restriction of phase space to regions of large q 2 also restores the validity of the OPE |7J and suppresses 
effects due to the details of the wave function of the b quark in the B meson. Taking only the region kinematically 
forbidden to b — > clP, q 2 > (tub — mo) 2 unfortunately introduces a low mass scale [H] into the OPE and the 
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uncertainties blow up to be of order {Aqcd /mi). Another drawback of this method is the elimination of higher 
energy hadronic final states, which may exacerbate duality concerns. 



1.1 Fermi motion and shape function 

Let us begin by discussing Fermi motion. This phenomenon, originally discovered in nuclear physics, is classically 
described as a small oscillatory motion of the heavy quark inside the hadron, due to the interaction with the 
valence quark; in the quantum theory it is also the virtuality of the heavy flavor that matters. Generally, as 
the mass of the heavy flavor becomes large, we expect that the heavy particle decouples from the light degrees 
of freedom and becomes "frozen" with respect to strong interactions. That is the basic assumption of HQET. 
However, even if Fermi motion can be neglected in the "bulk" of the phase space of the decay products, it 
still plays a role close to kinematical boundaries, such as the region of interest. Kinematically, that is easy to 
show, since a small virtuality of the heavy flavor in the initial state produces relatively large variations of the 
fragmentation mass in the final state. A b quark in a B meson has momentum 

V l = m b v" + (3) 

where w M is the four-velocity of the quark, which we can take at rest without any loss of generality: = 
(1; 0, 0, 0). fc M = p 1 ^ — m b v^ is the residual b quark momentum after the "mechanical" portion of momentum is 
subtracted off and it is of order Aqcd- If the momentum transfer to the final state leptons is q, the momentum 
and the invariant mass of the final state hadrons are 

p x = m b v^ + - q* 1 m x = p x (4) 

The boundary kinematic region is characterized by relations l|2|l: being /c M of order Aqcd, large values of Ex 
can originate only from m b v^ — q^ , inducing an almost light -like behavior 

mbv^-q* = (Ex,0,Q,E x ) + O(A Q cd) 

{m b v~q) 2 - 0{E x A QCD ). (5) 

The invariant mass of the final state hadrons is 

m x = (m b v + k - q) 2 = (m b v - q) 2 + 2k ■ (m b v - q) + 0(A QCD ) 

~ {m b v - q) 2 + 2E X k+ (6) 

where n M = (1, 0, 0, 1) and k + = k ■ n. Over most of phase space, the second term is suppressed relative to the 
first by one power of AQCD/m b , and so it may be treated as a perturbation. This corresponds to the usual 
OPE. However, in the region of interest the first two terms are of the same order. 
This can be also seen in a more compact way, imposing directly 

- q" = Exn" + fc"*, (7) 

where k^' is of order Aqcd- Then Eq. © can be written as 

m 2 x — (m b v + k - q) 2 = (m b v — q) 2 + 2k ■ (m b v — q) + 0{A 2 q CD ) 

~ E x k' -n + 2E x k-n = E x {k' + + k+). (8) 

A fluctuation in the heavy quark momentum of order Aqcd in the initial state produces a variation of the 
final invariant mass of the hard subprocess of order 

5m x ~0(A QCD Ex). (9) 

An amplification by a factor Ex has occurred, as anticipated. 

The differential rate in this region is therefore sensitive to the wave function f(k + ) which describes the 
distribution of the light cone component of the residual momentum of the b quark. The shape function is a 
non-perturbative function and cannot be calculated analytically, so the rate in that region is model dependent 
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even at leading order in Kqcd / m b- If we consider a heavy quark with the given off-shell momentum and a final 
state consisting of a massless on-shell quark at the tree level, we find for the shape function 

f(k + f art =5{k + + {m ^ Ex q? ) + O (a s ) , (10) 

as it should be, since m 2 x = 0. Selecting the hadronic final state, i.e. fc + , we select the light-cone virtuality of 
the heavy flavor which participates in the decay. 

We note that even with the amplification effect, Fermi motion effects are irrelevant in most of the phase 
space, where typical values for the final hadron mass are 

m\ ~ 0(E%), (11) 

in agreement with physical intuition. 

The shape function /(fc+) is written as a function of fc+ = fc° + ku; spatial components fcii and k± have 
been defined with respect to the mtv^ — direction (roughly the recoiling u quark). At this order, possible 
contributions due to k± are ignored. The shape function can also be seen as a resummation of the OPE to all 
orders in E x K QCD /m 2 x HOHEI]. 

Because the shape function depends only on parameters of the B meson, it is reasonable to expect that this 
leading order description holds for any B decay to a light quark. In general, we aspect QCD to factorize long 
distance effects into structure functions, with universal characteristics. In particular, one would like to get an 
estimation of structure functions for the semileptonic decays through the differential rate of B — > X s -f decay. 
However, due to the different kinematics between two body and three body processes, and to the presence of 
more than one energy scale (Ex, rrix and mj), universality of structure functions is not trivially applicable 

1.2 Perturbative resumming of large logarithms 

In general, the differential partial width is given by the convolution of a non-perturbative structure function with 
the perturbative calculable parton level differential distribution. Large remnants of the long distance dynamics 
occur also at the perturbative level by the presence of large logarithms near the threshold regions of phase space. 
Threshold resummation is a well known calculation technique which organizes the logarithmic enhancements 
to all orders in perturbation theory, thereby extending the QCD predictive power. In perturbative QCD, the 
hadronic subprocess in B — > X u lv consists of a heavy quark decaying into a light quark which evolves later into 
a jet of soft and collinear partons. The series of large infrared logarithms is due to an only partial cancellation 
of infrared divergencies in real and virtual diagrams. Let us consider f.i. the light quark produced in a process 
with a hard scale Q, evolving into a jet whose invariant mass is kinematically limited to a value m well below 
Q: the smaller integration region of the real diagrams induces a left-over double logarithm in the ratio Q/m. 
Multiple gluon emission occurs at high orders of perturbation theory, originating a double logarithmic series 



2 Radiative B decays 

Let us consider the radiative decay with a real photon in the final state, 

B -> X. + 7 (12) 

A systematic resummation is best done in the TV-moment space or Mellin space, because it leads to the expo- 
nentiation of the logarithmic corrections. In the ./V-moment space the threshold region corresponds to the limit 
N — > oo mil . Let us consider the Mellin transform of the normalized spectrum 

T R , N = f\l^t s ) N -^^dt s (13) 

Jo t_r dt s 
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In Mellin space threshold logarithms manifest themselves as a series in logiV. The order by order perturbative 
expansion contains double logarithmic contributions and it has the following schematic form 



Fr^n = 1+ a (c i2 L 2 + en L + cio) + a 2 (c 24 L A + c 23 L 3 + c 22 L 2 + c 21 L + c 20 ) H 

where L = logiV and a is the running coupling evaluated at the hard scale mj. 

The logarithmic terms have an exponential structure and we can write the following factorized form 

r_R,Ar = C R (a) a N (a) + d R , N (a), 

where the form factor <7/v in TV-space contains all the large logarithms to all orders in a 

oo n+1 

ox]> ^2^2g„ 



(14) 



(15) 



crjv(a) 



G N {a) 



„n t k 

a L 



n=l k=l 

= exp [Gi 2 a L 2 + Gna L + G 23 a 2 L 3 + G 22 a 2 L 2 + G 2 ia 2 L + G 34 a 3 L 4 + 
The coefficient function C R (a) has a perturbative expansion of the form 

C R (a) = 1 + aC R 1] + a 2 C ( R 2) + 0(a 3 ), 



(16) 



(17) 



(i) 

where C R are numerical coefficients, independent on N . The remainder function has a similar perturbative 
expansion; it depends on N, but it tends to zero in the threshold limit 



dn,N(ot) — * for N — > oo. 
The double sum in the exponent 116(1 is usually organized as a series of functions c/i 



a N (a) = exp L gx (/3 aL) + ^ a" g n+2 (foaL) 
The functions gi (A) have a power-series expansion: 

oo 

9i W = Yj9 mXn 



exp [Lg t (A) + 32 (A) + ag 3 (\) + a 2 g 4 (X) 



(18) 



(19) 



(20) 



where A = P^aL and [3q is the first coefficient of the (3 function expansion in a. By mantaining in the 
exponent of cr/v only the function g\ one approximates at leading order (LO); similarly keeping g\ and g 2 one 
approximates at next-to- leading order (NLO), g±, g 2 and 53 at next-to-next-to leading order (NNLO) and so on. 
This is equivalent to resumming series up to L(aL) n , (aL) n and a(aL) n , and so on. The explicit expressions 
of g, are known up to NNLO [13 [TBI EH EH 

The exponent in 1)16(1 can be obtained by means of the following resummation formula (141 1151 12()( 



G N (a) 



V N-1 



dz- 



1 - z 



2{1 - z) dk 2 1 

a -±A [a(k 2 )} + B [a(Q 2 (l - z))} + D [a(Q 2 (l - z) 2 )} \, (21) 



Q 2 (l-z) 2 



where Q = 2Ex is the hard scale of the process: in radiative decays, Q = mi,. The functions A (a), B (a) and 
D (a) have a standard fixed order expansion in a, with numerical coefficients 



A (a) = Aia + A 2 a 2 



B (a) = B ia + B 2 a 2 



, D (a) = Dia + D 2 a 2 



(22) 



A (a) describes the emission of partons which are both soft and collinear, B (a) describes hard and collinear 
partons, while D (a) describe partons which are soft and at large angles. The values of the known first coefficients 
of the functions A (a), B (a) and D (a) are reported in [TTj. 

By truncating a expansions for the functions A(a), B(a) and D(a) in eq. 1(21(1 one is implicitly assuming 
a <C 1. That is not always correct since the running coupling a (fcf) is integrated over all gluon transverse 
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momenta k t from the hard scale Q down to zero: inside the integration region, the Landau pole is hit and the 
running coupling diverges. A prescription has to be assigned to give a meaning to the formal expression (|21l) . 
Even after integration, an effect of the presence of the Landau pole persists: the series in eq. I|19|l is divergent as 
the higher order functions have factorially growing coefficients ^5J|2]. By truncating it to its first few terms, 
we stay within the so called fixed logarithmic accuracy. Moreover, the functions <?i(A) have branch cuts starting 
at A = 1/2 and going up to infinity. When A — > \ , by definition of A, a singularity in TV-space occurs at 
A — > exp [l/2f3oas (Q 2 ))], that is at N cr i t ~ Q/A (A is the QCD scale) 2 . The form factors are then formally 
well defined up to a critical value N C rit, above which they acquire a (completely unphysical) imaginary part. 
To return from A space to the physical space one uses the Mellin inverse transform: 



The evaluation of the inverse transform requires a prescription, in order to overcome the problems just mentioned 
and to obtain a form factor formally well defined in the whole t s -space, which resums all the logarithmic terms 
at the requested order. Several strategies can be pursued towards this result. The simplest possibility is 
to restrict oneself to a fiducial region in A-space below N crit . Another possibility is to use the so-called 
minimal-prescription |22| . which regularizes the form factor by means of an additional prescription for the 
contour integration in A-space in (|23|) . The problem of the presence of an integration over the Landau pole 
in resummation formulas, and of the ambiguities associated, has been also examined in the contest of the 
occurrence of infrared renormalons, to get information about non-perturbative effects, in the form of power- 
suppressed corrections to the cross sections (see f.i. HJ). 

Another approach is to give a prescription for the infrared pole directly in A-space, in such a way that 
the form factors arc well-defined for any A. Then it is not necessary to give a prescription for the Mellin 
inverse transform. A recent analysis |23| substitutes an effective coupling to the standard running coupling 
in resumming formula (|21|l . The effective coupling is built in a way to mantain the high energy behaviour 
of the standard running coupling, without reaching the Landau pole at low energies. Therefore the effective 
coupling never straddles far outside the perturbative domain. The resummation formula is free of Landau pole 
pathologies, and no prescription is needed because erjv(a) is analytic on the integration contour iffify . 

The effective form factor may also include absorptive effects related to the coupling constant in order to 
derive an improved expression for the resummation formula. Such absorptive effects are related to the decay 
of time-like gluons in the jet evolution. As well known from perturbation theory, at higher orders, one has to 
consider multiple emissions off the heavy and the light quark in B and secondary emissions off the radiated 
gluons. Because of the presence of these higher-order terms, the coupling in the resumming formula is evaluated 
at the transverse momentum of the primary emitted gluon |24|. By including the — iit terms in the integral 
over the discontinuity, i.e. the absorptive effects, usually neglected, the coupling in the resumming formulas is 
replaced by an effective coupling, evaluated at the transverse momentum of the primary emitted gluon |23II25| 



There is another advantage by using a perturbative approach with an effective coupling. In the standard 
approach, after resumming the large perturbative logarithmic contributions, one has to postulate a physically 
motivated non-perturbative model. Generally, an ad-hoc non perturbative form factor is convoluted with the 
perturbative distributions. Universal aspects of QCD dynamics, common to different processes, are not easily 
discovered this way. On the other hand, by describing different processes with the appropriate perturbative 
formulas, and the same effective coupling by assumption, one can deal simultaneously with perturbative and 
non-perturbative effects (without introducing model-depending parameters), avoid mistakes related to double 
counting, and underline universal effects by comparing with the data |23l 1251 12(i| . 

2 Let us also observe that the degree of singularity of the functions gi for A — * 1/2, and therefore also of the form factor, increases 
with the order of the function, i.e. with i 1211 . 




(23) 



a -> S(fc 2 ). 



(24) 



5 



3 Universal aspects of QCD dynamics 



Let us now consider the decay 

B -> X u + I + v. (25) 

It is possible to obtain a factorized form for the triple differential distribution, the most general distribution in 
process <|2l)|) (its integration leads to all other spectra): 



1 d 3 T 



r dxdudw 



C [x, w\ a(w m,b)] & [u; a(w nib)] + d[x, u, w,a(wnib)] , (26) 



where: 



and 



w = ^E. (o < w < 2), x = — (0<x<l) (27) 



E x - yj& x - m 2 x _ l-y/l-(2mx/QY ^ /mx 
E X + VE X ~ m\ 1 + L _ {2mx/Q f " V Q 



(0<u<l). (28) 

We have called Q the hard scale of the process and, at the threshold, set Q = 2 Ex ■ Both the logarithms 
and the argument of the running coupling depend on the kinematics of the problem, by means of the hard 
scale Q. There are two important kinematical differences with respect to the radiative case. First, in the three 
body semileptonic decay, the distribution also depends on the charged lepton energy E\. Second, while in the 
radiative decay it is always Q — 2 Ex = nib, m process (I25f) there are regions of phase space where Ex is 
substantially reduced 3 . 

In analogy with what done in the previous section, we can study the distribution in the Mellin space, defining 

a N (a) ee / duil-u^- 1 a(u- a) (29) 



At this level, there is universality among radiative and semi-leptonic decays, meaning that the same QCD form 
factor on appears in both distributions Ijl5|l and l|29(l . Consequently, the form factor in the physical space is 
the same cr(u; a). It is evaluated at the argument u ~ m x /(AE x ) in the semileptonic case. In the radiative 
case, by imposing the kinematical relation between hadronic energy Ex 3 and hadronic mass mx s , we have 



\Ex s =m b /2(l+m 2 x /mf) 



t s . (30) 



where t s = m x /m 2 . 

The coupling constant argument is set at the hard scale Q = 2E X in both processes; in the radiative decay, 
that implies it is fixed to mb. 

This simple connection between radiative and semileptonic processes is sometimes lost when one passes 
to the double or single differential distributions for the processes 112 5 ft ^1^2 El ED- All double and single 
distributions are obtained by integrating the triple differential distribution (|2*fJ|l . As seen before, in semileptonic 
decays the hadronic energy E x is not fixed and it can be integrated over; this affects the infrared structure of 
the distribution, since the form factor a depends on E x |181 1111 IT21 This class of distributions, f.i. the 
single differential distribution in m x , cannot be directly compared with the radiative spectra. The structure of 
the threshold logarithms is not the same and there is no universality of long distance effects. On the other side, 
distributions as the single differential one in E x , where the energy E x is not integrated over, keep the infrared 
structure of the radiative spectrum. They can directly be related via short-distance coefficients to the radiative 
spectrum. 



3 One can think, f.i., to the kinematical configuration with a large invariant mass for the lepton pair 
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